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C/} ! We describe the crystalline phase of color superconducting quark matter. This phase 
\Q ', may occur in quark matter at densities relevant for compact star physics, with possible 
implications for glitch phenomena in pulsars [ |J]. We use a Ginzburg-Landau approach 
7—i • to determine that the crystal has a face-centered-cubic (FCC) structure [ 0] . Moreover, 
our results indicate that the phase is robust, with gaps, critical temperature, and free 
energy comparable to those of the color-flavor-locked (CFL) phase Our calculations 
also predict "crystalline superfluidity" in ultracold gases of fermionic atoms [[|. 

Cold dense quark matter is a color superconductor [||. At asymptotically high densities, 
the ground state of QCD with quarks of three flavors (it, d, and s) is the color-flavor- 
locked (CFL) phase [||. This phase features a BCS condensate of Cooper pairs of quarks 
that includes ud, us, and ds pairs. At intermediate densities, however, the CFL phase 
can be disrupted by any flavor asymmetry (such as a chemical potential difference or 
^| a mass difference) that would, in the absence of pairing, separate the Fermi surfaces. 
In the absence of pairing, electrically neutral bulk quark matter with m U ( i = and 
m s 7^ features a nonzero electron density (/i e ~ m^/4/i) and three disparate quark 
Fermi momenta: p F ~p F + m^/4/i, p s F « p u F — m 2 s /AyL. (Note that decreasing /i enhances 
the flavor disparity.) Accounting for pairing effects modifies this picture: starting in 
the CFL phase at large fi, as we decrease \i the CFL phase remains "rigid" [ ||, with 
coincident quark Fermi surfaces and no electrons, until either hadronization or a first- 
order unlocking transition, whichever comes first. Unlocking occurs at fi w m^/4A , 
where Ao is the CFL gap. Its value and that of m s are density dependent and sufficiently 
uncertain that we do not know whether unlocking occurs before hadronization. Here, we 
pursue the consequences of assuming that unlocking occurs first. 

In quark matter below the unlocking transition, pairing can still occur. One option 
is single-flavor pairing (uu, dd, ss), but these J = 1 condensates have very small gaps [ 
|[. Crystalline color superconductivity is more robust [ [l], §]. We propose that unlocked 
quark matter is in the crystalline phase, which therefore occupies the window in the QCD 
phase diagram between the CFL and hadronic phases (Fig. 1). The crystalline phase 
has the unique virtue of allowing pairing between quarks with unequal Fermi surfaces. 
It was originally described by Larkin, Ovchinnikov, Fulde, and Ferrell (LOFF) [| as a 
novel pairing mechanism for an electron superconductor with a Zeeman splitting between 
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Figure 1. Schematic QCD phase diagram 
with quark-gluon plasma, crystalline (XTAL) 
and CFL color superconducting phases, the lat- 
ter separated by the unlocking transition. 



Figure 2. A unit cell of the FCC crys- 
tal, with contours where the gap func- 
tion is positive (black), zero (gray), and 
negative (white). 



spin-up and spin-down Fermi surfaces, neglecting all orbital effects of the magnetic per- 
turbation. Quark matter is a more natural setting for the LOFF phase, as it features a 
"flavor Zeeman effect" with no orbital complications. Cooper pairs in the LOFF phase 
have nonzero total momentum: a quark with momentum p is paired with a quark with 
momentum — p + 2q. The magnitude |q| = qo is determined by the separation between 
Fermi surfaces while the direction(s) q is (are) chosen spontaneously. Each quark in a pair 
resides near its respective Fermi surface, so pairs can form at low cost and a condensate 
results. For a given q, the quarks that pair are only those in "pairing rings", one on each 
Fermi surface; these circular bands are antipodal to each other and perpendicular to q. 

If each Cooper pair in the condensate carries the same total momentum 2q, then in 
position space the condensate varies like a plane wave: (^(x)'0(x)) ~ Aexp(2iq • x) 
meaning that translational and rotational symmetry are broken. If the system is unstable 
to the formation of a single plane-wave condensate, then we expect that a superposition 
of multiple plane waves is more favorable: (?/>(x)?/>(x)) ~ i q i =?0 A q e 2iqx . Each A q 
correponds to condensation of Cooper pairs with momentum 2q, i.e. another pairing ring 
on each Fermi surface. As we add more plane waves, we utilize more of the Fermi surface 
for pairing, with a corresponding gain in condensation energy. On the other hand, the 
rings can "interact" with each other: condensation in one mode can enhance or deter 
condensation in another mode. We need to compare various different choices for the set 
of q's whose A q 's are nonzero; that is, we need to compare different crystal structures. 
To this end, we calculate the Ginzburg-Landau free energy f2({A q }) for the LOFF phase. 

We consider a simplified model with massless u and d quarks only, and introduce chem- 
ical potentials \i u = p> — 5fi, \id = fi + Sfi. In this toy model, we vary 5/i by hand to emulate 
the effect of splitting the Fermi surfaces as described above. (The LOFF phenomenon 
can equally well be induced by mass differences between quarks [ [K| . ) We work at zero 
temperature, as appropriate for compact star physics. We model the interaction between 
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quarks with a pointlike four-fermion interaction that mimics single gluon exchange. 

At 5fi = 0, the system forms a BCS superconductor with a gap A . As we increase <5/i, 
the system exhibits a "rigidity" analogous to that of the CFL phase: despite the imposed 
stress 5fi, the gap stays constant and the Fermi surfaces remain coincident. When 5fi 
reaches a critical value 5^ = A /\^2 « 0.707A , the BCS phase "breaks" and the Fermi 
surfaces separate (this is analogous to CFL "unlocking" ) . For 5ji > 5fi 1 a crystalline phase 
is possible. If we consider the simplest LOFF ansatz, a single plane wave condensate [ |TJ, 
we find that the LOFF state is favored in an interval 5fi 1 < Sfi < <5/i 2 ; at Sfi l there is a 
first-order transition from BCS to LOFF, while at 5/i 2 ~ 0.754A there is a second-order 
transition from LOFF to the normal state (unpaired quarks). Each Cooper pair has a 
total momentum 2q, where |q| = q « 1.205/i; this implies that the pairing ring on each 
Fermi surface has an opening angle ipo ~ 2 cos _1 (<5/z/go) ~ 67.1°. 

Motivated by the second-order transition at 5fi 2 , we write a Ginzburg-Landau free 
energy valid in the small order parameter limit [0. The most general expression consistent 
with translational and rotational symmetry is 

fioc «A;A q +^J(n)A; i A q2 A; 3 A q4 +i^K(o)A; i A q2 A; 3 A q4 A; 5 A q6 +---(l) 

q,|q|=<JO □ O 

The symbol □ represents a set of four equal- length vectors (qi, q 2 , q3, q4), |q«| = qo, with 
<li — c[2 + ci3 — — 0, forming a closed figure. Similarly, the symbol O represents a set of 
six equal- length vectors (qi, . . . , q 6 ) with qi — q2 + q3 — q4 — qs + q6 = 0. The quadratic 
coefficient a changes sign at 5/i 2 showing the onset of the LOFF plane-wave instability: 
a ~ (6fi — <5/z 2 ) / '6^2- Analysis of the quadratic term shows that all modes on the sphere 
|q| = go become unstable at 5fi < <5/i 2 , and that at larger 5/i it is these modes that are 
the least stable. The quartic and sextic coefficients J and K characterize the interactions 
between modes, and thus determine the crystal structure. For candidate crystal structures 
with all A q 's equal in magnitude, we can evaluate aggregrate Ginzburg-Landau quartic 
and sextic coefficients j3 and 7 as sums over all rhombic and hexagonal combinations of 
the q's: (3 = Y^,u 7 = Xlo ^(P)- Then for a crystal with P plane waves we obtain 

Q(A) oc PaA 2 + \l3A A + I7A 6 + • • • and we compare crystals by calculating (3 and 7 to 
find the structure with the lowest Q. 

We have analyzed many candidate crystal structures in [ 0; here, we summarize the 
lessons learned. Recall that each q corresponds to a pairing ring on each Fermi surface 
with opening angle ipo ~ 67. 1°. The first lesson is that adding more pairing rings lowers the 
free energy only as long as rings do not intersect. No more than nine rings can be arranged 
on the Fermi surface without any overlaps. The second lesson is that "regular" structures 
(those with nonintersecting rings whose wave vectors can be assembled into closed figures 
in many ways) are combinatorically favored. There are no particularly regular nine-wave 
structures, but there is a very regular eight-wave structure: eight wave vectors q pointing 
towards the eight corners of a cube, forming the eight shortest vectors in the reciprocal 
lattice of a face-centered-cubic (FCC) crystal. This FCC crystal has by far the lowest 
Ginzburg-Landau free energy of all the structures we considered. The gap function is 
(V>(x)V>(x)) ~ A[cos^(x + y + z) + cos^(x-y + z)+cos^(x + y-z)+cos^(-x + y-z)] 
with lattice constant a = V3ir/\q\ m 6.012/A . A unit cell is shown in Fig. 2. The effective 



field theory for the phonons of this crystal has recently been constructed [ [TT| . 
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The Ginzburg-Landau free energy for the FCC structure has coefficients (3 and 7 that 
are large and negative. This guarantees that the phase transition to the crystalline phase 
cannot be a second-order transition at 5ji = 5{i 2 and must instead be a strongly first-order 
transition at some Sfi = 8fi^ 3> 5/i 2 , at which the magnitude of the gap jumps from zero to 
a large nonzero value. Were we to push beyond sextic order to obtain a bounded Ginzburg- 
Landau free energy, the resulting calculation of the gap and ground state free energy 
would be quantitatively uncontrolled because the order parameter is nowhere nonzero 
and small. Qualitatively, the large negative coefficients indicate a robust crystalline phase, 
with gaps and free energy comparable to those of the BCS phase. Furthermore, because 
the transition to the crystalline phase occurs at 5fi^ ^> Sfj, 2 , the crystalline phase is favored 
over a wide 6fi interval. Our results thus indicate that wherever unlocked quark matter 
occurs, it will be in the crystalline phase with an FCC crystal structure. 

The power of the Ginzburg-Landau formulation is that it organizes the calculation in 
such a way that simple qualitative lessons about what features make a crystal structure 
favorable emerge. We immediately learn that, among all the possibilities, the FCC struc- 
ture of Fig. 2 is the best choice. This is the first time (in either the condensed matter 
or the QCD incarnations of this subject) that sufficiently general three dimensional crys- 
tal structures have been compared. The structure that we find to be most favorable, 
although simple with hindsight, was never considered in the classic condensed matter 
physics papers. Unfortunately, but not unusually, the Ginzburg-Landau calculation pre- 
dicts a first-order phase transition and in so doing predicts its own quantitative demise. 
With a qualitative understanding of what crystal structure wins (that is, with confidence 
in what ansatz to make) we are currently doing the variational calculation needed to de- 
termine the gap, free energy, and 5fi^ quantitatively within our toy model. Going beyond 
our toy model, a full three-flavor analysis is desirable, with ud, us and ds crystalline con- 
densates and a strange quark mass. In an atomic physics context, the FCC condensate 
may be detectable as a "crystalline superfluid" in an ultracold gas of fermionic atoms [ [|. 
In an astrophysics context, with the crystal structure known a calculation of the pinning 
force can now proceed, in order to determine whether (some) observed pulsar glitches 
originate from pinning of rotational vortices at intersections of crystal nodal planes. 
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